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Abstract 

We consider the simplest nontrivial super symmetric quantum mechanical system 
involving higher derivatives. We unravel the existence of additional bosonic and 
fermionic integrals of motion forming a nontrivial algebra. This allows one to obtain 
the exact solution both in the classical and quantum cases. The supercharges Q, Q 
are not Hermitially conjugate to each other anymore, which allows for the presence of 
negative energies in the spectrum. We show that the spectrum of the Hamiltonian is 
unbounded from below. It is discrete and infinitely degenerate in the free oscillator- 
like case and becomes continuous running from — co to cxd when interactions are 
added. Notwithstanding the absence of the ground state, there is no collapse, which 
suggests that a unitary evolution operator may be defined. 

1 Introduction 

It was suggested in Refs. [1,2] that the Theory of Everything may represent a conven- 
tional supersymmetric field theory involving higher derivatives and living in flat higher- 
dimensional space. Our Universe is associated then with a 3-brane classical solution in 
this theory (a kind of soap bubble embedded in the flat higher-dimensional bulk), while 
gravity has the status of effective theory in the brane world-volume. 

Generically, higher-derivative theories involve ghosts [3] described usually as negative 
residues of the propagator poles and/or indefinite metric of Hilbert space. Speaking in 
more direct physical terms, the presence of ghosts means the absence of the lower bound 
(the ground state) in the spectrum of the Hamiltonian. This more often than not leads 
to violation of causality or unitarity or both (see e.g. the recent discussion in [4]). 

The problem of ghosts was discussed recently in Refs. [1,5,6]. In particular, in Ref. 
[5] a nontrivial quantum mechanical higher-derivative system was presented where the 

*0n leave of absence from ITEP, Moscow, Russia. 
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spectrum was bounded from below and hence the ghosts were absent. To be more precise, 
the spectrum of this system has no bottom in the free "Pais-Uhlenbeck oscillator" case, 
but the bottom appears as soon as the interaction (of a certain kind) is switched on. 
When the interaction constant a is small, the ground state energy behaves as —C/a. 
Negative and large by absolute value, but finite. 

This example was not supersymmetric, however, and the mechanism by which the 
ghosts were killed there seems to be specific for nonsupersymmetric systems. In Ref. [6], 
we considered a supersymmetric model {bD superconformal gauge theory reduced to 0+1 
dimensions) which naively involves ghosts. But we showed that one can effectively get rid 
of them, if working in reduced Hilbert space where the Hamiltonian is Hermitian and its 
spectrum is bounded from below. It is supersymmetry which helps one to do it. Indeed, 
the standard minimal supersymmetric algebra 

= = 0, 

{Q,Q]^2H, (1) 

where the supercharges g, Q are Hermitially conjugate to each other, implies that all 
eigenvalues of the Hamiltonian are non-negative and the ground state with zero or positive 
energy exists. 

Though ghost-ridden, the model of Ref. [6] did not involve higher derivatives in the 
Lagrangian. The motivation of the present study was to find out whether the ghost-killing 
mechanism found in [6] works also for higher-derivative supersymmetric theories. To this 
end, we considered the simplest higher-derivative supersymmetric quantum mechanical 
system with the action 



dtdOdO 



'-{VX)j^{VX) + V{X) 



(2) 



{X is a real supervariable) . We will be mainly interested in the case where V{X) is a 
polynomial having the form 

V(X) = -^-^. (3) 

We found that, though certain technical similarities between this system and the system 
considered in [6] exist, the physics in this case is essentially different. In particular, 
there are no compelling reasons to censor the negative energy states out of the spectrum. 
However, in spite of their presensc (so that the spectrum is unbounded both from above 
and from below), this does not lead to disaster. Irreversable loss of unitarity is usually 
related to collapse phenomenon where singularity is reached in finite time. In our case, 
there is no such collapse. Moreover, the eigenstates of the Hamiltonian that we found 
have real energies and it is reasonable to expect that the unitary evolution operator can 
be defined. 

In the next section, we describe the model, write down the component expressions for 
the Lagrangian, supercharges and the Hamiltonian. We discuss the trivial noninteractive 
case V{X) oc X"^ and then the generic case. We exhibit the presense of certain additional 
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integrals of motion, which makes the problem exactly soluble. In Sect. 3, we discuss 
the classical dynamics and, for the superpotential ([3]), write the solutions to the classical 
equations of motion explicitly. We show that there is no collapse and the solution exists 
at all times. It has an oscillatory behavior with linearly rising amplitude. In Sect. 4, 
we address the quantum problem and find the spectrum and the eigenstates. In Sect. 5, 
we consider a more complicated system where the higher-derivative term is added to the 
conventional kinetic term. Its classical dynamics is even more benign than the dynamics 
of pure higher- derivative theory — the amplitudes do not rise linearly anymore and the 
motion is bounded in a finite region of the phase space. The spectrum of the mixed system 
is probably discrete, but dense everywhere. In Sect. 6, we discuss briefly a model where 
still extra time derivative is added. The last section is devoted as usual to concluding 
remarks and speculations. 



2 The model. 

Let us express the action ([2]) in components. To this end, we substitute there 

X = X + eip + ipe + Dee , 

^ 9 .T^d - d .^d 
V = — + ie—, V = =-«6'— , 

de dt' de dt ' 

and integrate over dede. We obtain 

L = xD + V'{x)D + V"{x)iljilj + iljij . (4) 

Note that this Lagrangian involves twice as much physical degrees of freedom compared 
to standard Witten's supersymmetric quantum mechanics [7], 



^ stand 



dede 



^vxvx + V{X) 



x^ + D^ 



+ tpij'^ + V'{x)D + V"{x)^^ . (5) 



Indeed, the field D enters the Lagrangian (jl]) with a derivative and becomes dynamical. 
In addition, does not coincide anymore with the canonical momentum of the variable 
but represents a completely independent complex fermion variable not necessarily 
conjugate to ip- It is convenient to denote it by x s-nd reserve the notation ip, x for the 
canonical momenta 

X = 'i'Px = ^p = ip^ = -ix . 

Introducing also 

p = p^ = t); P = Pd = x , 

we can derive the canonical Hamiltonian 

H =pP- DV'{x) +ipx- V"{x)xiJ . (6) 
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The Lagrangian (jlj) (with x substituted for ip) is invariant (up to a total derivative) with 
respect to the supersymmetry transformations, 



5eX = ex + ipe , 
S^ifj = e{D — ix) , 
SeX = KD + ix) , 

6,D = z(ex - i^e) . (7) 
The corresponding Nother supercharges are 

Q = ^[p + ^V'{x)]-x{P-^D) , 

g = ^(P + iD) - x[p - tV'ix)] . (8) 

One can be convinced that the algebra ([1]) holds, but, in contrast to the standard SQM, 
Q and Q are not Hermitially conjugate to each other. This is the main reason for all the 
following complications. 

Consider the simplest case, 

V{X) = . (9) 

It is convenient to make a canonical transformation 

x+ + X- I — — , , I — — , , p+ — p_ 



i; = t±^^ ^ = t_^^ ^ = V^(V;+ + V'-), x=v^(^--^+). (10) 



In terms of the new variables x±^p±,il^±,'4)±, the supercharges and Hamiltonian acquire 
a simple transparent form 

Q = ip+{p+ - iujx+) + ip-{p^ - iujx^) = + Q_ , 

Q = 'ip+{p+ + iujx+) - ip^ip^ + iujx^) = - ; (11) 

g = ^ ^ + +^^+^+- ^- ^ - . (12) 

In other words, the system represents a combination of two independent supersymmet- 
ric oscillators such that the energies of the second oscillator are counted with the negative 
sign (In nonsupersymmetric case, such a system was first considered by Pauli back in 
1943 [8]). The states are characterized by quantum numbers {n±,F±}, where n± are 
nonnegative energies characterizing the excitation levels of each oscillator and F± = 0, 1 
are the fermion numbers, the eigenvalues of the operators The spectrum of the 

Hamiltonian 

En+,n_ = a;(n+ - n_) (13) 
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is infinitely degenerate at each level depending neither on n+ + n_ nor on F±. The 
spectrum (fT3l) is discrete involving both positive and negative energies. 

We see that, in spite of supersymmetry, the spectrum has no bottom and hence involves 
ghosts. In contrast to what was the case for 5D superconformal theories [6] , the negative 
energy states have the same multiplet structure as the positive energy ones and there 
are no "scientific" reasons (i.e. the reasons based on certain symmetry considerations) to 
exclude these states from the spectrum. 

However, these ghosts are definitely of benign variety. Actually, when the system 
consists of several noninteracting subsystems whose energies are individually conserved, 
the sign with which these energies are counted in the total energy is a pure convention. 
The problems may (and do usually) arise when the subsystems start to interact. Then, 
if it is the difference rather than the sum of the energies of individual subsystems that 
is conserved, there is a risk that the individual energies would rise indefinitely leading to 
the collapse with associated unitarity and causality loss. 

What happens in our case ? A proper way to include interactions is to modify the 
superpotential The key observation is that for any superpotential V{X) the system 
involves besides H,Q,Q two extra even and two extra odd conserved charges. They can 
be chosen in the form 

p2 

N = —- V{x) , 

F = ipip-xx , 
T = ij[p-iV'{x)]+x{P + iD) , 

f = ^{P -iD) + x[p + iV'{x)] . (14) 

The superalgebra {H, N, F; Q, Q, T, T) has the following nonvanishing commutators: 

{Q,Q} = {T,f} = 2H- 
[Q, F] = Q, [Q, F] = -Q, [T, F] = -T, [T, F] = T; 

[Q, N] = [T, N] = -IQ, N] = [f , N] = ^ . (15) 

Now, T and T are the extra supercharges, the subalgebra involving the operators {H-Q, Q, T, T) 
coincides with the standard subalgebra of extended J\f = 2 supersymmetry S2. □ This 
leads to 4-fold degeneracy of each nonvacuum level in quantum problem (but does not 
lead necessarily to positivity of their energies as Q is not conjugate to Q and T is not 
conjugate to T). F is the operator of fermion charge. As defined, it takes values for the 
states with the wave functions \1/ oc 1 and \1/ oc ipx, the value 1 for the states oc ip and 
the value —1 for the states \E' oc x- (The convention is somewhat unusual, but one could 
bring it to the standard form by interchanging x and x-) Finally, the conserved charge 
iV is a new object that is specific for the problem in hand. The corresponding symmetry 
of the action is 

D ^ D + ax . (16) 



is an ideal of the superalgebra (fT5)) and hence the latter is not simple. It represents a semidirect 
sum of the Abelian Lie algebra {F, N) and 52 . 
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Indeed, the Lagrangian (jlj) is shifted by a total derivative after this transformation. 



3 Classical dynamics. 

Let us disregard the fermion variables and concentrate on the dynamics of the bosonic 
Hamiltonian 



It involves two pairs of canonic variables. The presence of the extra integral of motion 
implies that the system is exactly soluble and seems to imply that the variables can be 
separated and the classical trajectories represent toric orbits. This would be true for a 
conventional system with positive definite kinetic term. In our case, the situation is more 
complicated. To begin with, the variables cannot be easily separated. Indeed, excluding 
the momenta from the corresponding canonical equations of motion, we obtain 



The equation for x does not depend on D, but the equation for D does depend on x for 
generic V{x). 

Let us try first to add the cubic term to the superpotential V{X). As we see, the 
same function taken with the negative sign plays the role of the potential for the variable 
X. If V{x) oc x^ at large x, the potential is not binding and the motion is infinite such 
that infinity is reached at a finite time. This is the collapse signalizing the presence of 
the ghost of malignant variety. 

Let us choose now the polynomial negative definite superpotential. It involves only the 
even powers of X. The simplest nontrivial case is Eq.(l3]). The potential is confining now 
and the equation of motion has a simple solution representing an elliptic cosine function 
with the parameters depending on the integral of motion A^, 




(17) 





(18) 



x{t) 



Xq cn[f2t, k] 



(19) 



with 




(20) 



Here k is the parameter of the Jacobi elliptic functions. [9] 



^ Recall that, if k gJO, 1[ and t — —===== ^ then the elliptic functions are: snt = sin0, 




cnt = COS0, dnt = \Jl — k'^ sin^ (j). The functions sn, cn, dn are periodic with period AK where 



" Vl-fe^ sinks' 
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The equation for D represents an eUiptic variety of the Mathieu equation. Generically, 
it is not the simplest kind of equations, but in our case the solutions can be obtained in 
a rather explicit form. [f| One of the solutions is 



Di{t) oc x{t) =oc sn [fit, k] dn [fit, k] . 



(21) 



The second solution can be found from the condition that the time derivative of the 
Wronskian W = D1D2 — -D2-D1 vanishes. We find 



D2{t) oc x{t) 



dt' 



=oc sn [fit, k] dn [fit, k] 



dt' 



x2(t') ™ '"J -J J sn2 [fit', k] dn^ [fit', k] 

When u; = 0, the integral in Eq. fl22l) can be done analytically and we obtain 

D{t) = Asn fit, ^172 dn fit, v^l72 + 
fit, ^172 dn fit, ^172 } 



(22) 



Bicn 



fit, Vl/2 



fit sn 



(23) 



Two independent solutions ( |2T]) . ( l22l) exhibit oscillatory behaviour with constant and 
linearly rising amplitude, correspondingly (see Fig.l). The energy E does not depend 
on A and is 



E = BX^/\4N) 



3/4 



(24) 




Figure 1: The solution of the equation (ITSl) for D{t) with the parameters u = X = N = 1 
and inital conditions -D(O) = 1,D'{0) = 0. 



^We thank N. Nekrasov for this remark. 

* One can remind the situation for the ordinary Mathieu equation. In generic case, its solutions, the 
Mathieu functions, exhibit osciUatory behaviour with the amphtude that either osciUates itself or rises 
exponentially. But for some special characteristic values of parameters, the amplitude stays constant or 
rises lineary, like in our case. 
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4 Quantum dynamics 



4.1 Bosonic system 

Consider first the bosonic Hamiltonian (fT7|l . Let us prove that the corresponding evolution 
operator is unitary. To this end, it is convenient to perform a partial Fourier transform 
and consider the wave function in the mixed representation, 

/oo 
e-'^^^(x,D)rfD, (25) 
-oo 

The Schrodinger equation 

at 

{Kb is the operator obtained from Hb by the corresponding canonical transformation) 
for the function \l/(x, P) (we will not right tildas anymore) represents a linear first order 
differential equation. 

This equation can be easily solved by the charateristics method %. The characteristic 
system is here 

X = P 

P = V'{x) (27) 

The equations (!27j) represent a half of original Hamilton equations of motion for the 
system (1171) . They can be interpreted as the Hamilton equations for the system described 
by the "Hamiltonian" — V{x). The latter coincides with the extra integral of motion 
N{P,x) defined before and should not be confused with the true Hamiltonian Hb- 

Let us denote by F* the flow determined by (1271) . By definition, we have F*(xo,Po) = 
{xt,Pt). We clearly see that the Schrodinger equation fl26|) is solved by 

^,(x,P)= ^o(r-*(x,P)) (28) 

with an arbitrary "^qIxjP). Moreover, as —V is confining, the flow F* is well defined 
everywhere in for all times and this property entails that the Hamiltonian Kb and 
hence Hb are essentially self-adjoint. 

At the next step, we will solve the stationary spectral problem for Hb and find the 
eigenstates. We will construct the states where not only the Hamiltonian Hb, but also 
the operator have definite eigenvalues. The system is integrable and a regular way to 



^see for example Ref. [11] 

^this means that Kb, Hb have a unique self-adjoint continuation in L^(R^) starting from the space 
of smooth, finite support functions on R^. 

note for purists: as most of these states belong to continuum spectrum, they represent generalized 
eigenstates of Hb and N. 
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solve it is to go over into action-angle variables. There is some specifics in our case. We 
will follow the standard procedure not for Hb (it is not possible as the variables cannot 
be separated there) but for the quasi- Hamiltonian N involving only one pair of variables 
(P, x). Thus, we perform a canonical transformation S: {x,P) ^ {I is the action 

variable, ip is the angle, / G ]0, +00 [, ip e [0, 2n[) such that in this new coordinates system 
the flow is 

T\S-\l,cp)) = S-\l,ip + tcriI)) (29) 
where a{I) = dN/dl. Let us recall that / is given by the following integral 

I = (p Pdx = 77- / dxdP , 

27r J 27r Jni^x,p)<No 

where A^'o is the energy coinciding in our case with the value of the integral N on the 
trajectory. For the potential ([3]), one can derive 

with f2, k written in Eq. (!20|) . In the purely quartic case, c<j = 0, A = 1, 



with ^ 

K = K{1/V2) = ^^-^^^1.85. 

The explicit expressions for the canonical transformation 5* from the action-angle 
variables to the variables x, P are in this case 

X = cn {p 

V 

P = -fi2(j)sn(^^<^^dn(^^^^. (32) 

with the angle ip G M/27rZ and the positive action / > 0. 

In the representation where the wave function depends on / and ip, the solution 
fl28l) to the Schrodinger equation takes the form 

^^(J, ip) = f/(t)^o(/, = *o(/, f - ta{I)) . 

In this representation, U(t) is a unitary evolution in the Hilbert space L'^(]0, oo[, R/27rZ). 
Its generator is a new quantum Hamiltonian: 

nij = -ta{I)—. (33) 
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The Hamiltonians Hb and Kb are unitary equivalent to the Hamiltonian Ti,. 

Using a Fourier decomposition in the variable (f, we have an explicit spectral decom- 
position for H. If ^(/, ^) = J2 ^n(/)e™'^, then 

Hi^il, ^) = J2 ^o^(^)^n(/)e*"^. (34) 
Substituting in = na{I) the expression fl30l) . we derive the quantization condition 

1/4 

(35) 

In the limit A — 0, the dependence of the left hand side of Eq.( l35|l on disappears and we 
reproduce the simple oscillator quantization condition E = ujn coinciding with Eq.( fT3l) . 
When A 7^ 0, the right hand side of Eq.(l35]) depends on N [En ~ T^n{\Nf/^/y2K{l/^)] 
for large A^] and only a certain combination of E and N is quantized, but not the energy by 
itself. For illustration, the function Ei{N) is plotted in Fig|2]for two choices of parameters. 
The dependence of cr on J and hence En on N reveals that the spectrum is continuous 



En 



irn 
2K{k) 




Figure 2: The dependence Ei{N). The lower curve corresponds to the choice u = 1, X = 
0.1 and the upper one to u = X = 1. 

here, with eigenvalues lying in two intervals ] — oo, — u;] U [a;,+oo[ plus the eigenvalue 
{0}. The same qualitative picture (continuum spectrum which can be supplemented 
by isolated eigenvalues) holds for generic binding potentials —V{x), in particular, for 
generalized anharmonic oscillators, V{x) = —aox"^^ + aix^^~^ + ■ ■ ■ + a2i, > 0, i > 1, 
where (t(/) ~ for large I. 
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The generalized eigenfunctions of the Hamiltonian (!33|) are labelled by the parameters 
Jo e M and n G Z, 



^i,n{I,f) = 5(/-/o)e*"^ . (36) 
Going back to the original variables using Eqs. fl25ll32p . we obtain 

^en{x,D) = ] e-^(-'^), (37) 

+ V{X) 

where 

S{x, D) = D^2[N + V{x)] + ^ r , ^ ^ (38) 

is nothing but a classical action function of the original system [not to confuse with the 
constant I proportional to the action on a closed trajectory of the reduced system ([2] 
S{x, D) satisfies a system of generalized Hamilton- Jacobi equations 

1 f dSV 



For the superpotential ([3]), the second term represents the elliptic integral of the first kind, 
E f = ISLf farosm ( ^] , -J^) . (40) 

with xo and k given above. It is convenient to express it into inverse elliptic cosine function 
arccn{u, k). Substracting an irrelevant constant, we may rewrite Eq. fl38|) as 



X 



S{x,D)=D^2N-uj^x^-^-^^^^=^aTccn\^-^,kj . (41) 

For E,N satisfying the quantization condition (j35l) . the wave function (!37j) is single- 
valued. 

We see that the exact solution (1371) differs from the semiclassical wave function e**^ by 



the extra factor 1/ \/2N — lu'^x'^ — Xx^/2. For large enough x and nonzero N, the function 
falls down exponentially (we have to choose the sign of the square root in accordance with 
the sign of D). For intermediate x, the function oscillates in D and behaves as a plane wave 
continuum spectrum solution. When V{x) = N, the wave function involves a singularity, 
with the normalization integral diverging logarithmically at this point. 
Two natural questions are in order now. 

1. Is this singularity at finite value of x dangerous ? 
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2. How come the non-normalizable wave functions (!37|) describe also the zero energy 
states ? The point E = is isolated and one expects that the eigenfunctions with 
zero energy belong to L^. 

Let us answer first the second question. There are infinitely many states of zero 
energy. In the action-angle variables, any function g{I) —>■ g{N) not depending on is an 
eigenfunction of (l33l) with zero eigenvalue. In original variables, this gives the function 

^o{x, D) = {2n)-'/' 1^ al^Y' ^^"^0 "'^'^'''^^ ■ ^^^^ 

The solution (!37|) is obtained, if substituting in Eq. p2!) ^(A^) = S{N — Nq). But we may 
also choose the basis gk{N) = N''e~^ , k = 0, 1, . . . (its orthogonalization gives the La- 
guerre polynomials) giving the normalized zero-energy solutions without any singularity. 
Any smooth function can be expanded into this basis. The distribution 6{N — Nq) can, 
of course, be represented as a limit of a sequence of smooth functions. 

The existence of the normalized zero energy states together with continuum states 
could somewhat remind the maximal supersymmetric Yang-Mills quantum mechanics [10]. 
There are two differences: (i) The latter is a conventional supersymmetric system and 
the zero-energy states have the meaning of the vacuum ground states; (ii) In our case for 
a; 7^ 0, the zero energy state is separated by a gap from continuum. For the maximal 
SYM quantum mechanics, there is no gap. 

The inverse square root singularity of the continuum spectrum functions \l/ e^o has the 
same nature as the divergence of their normalization integral at large D. It is benign and 
physically admissible. Indeed, the physical requirement for the systems with continuum 
spectrum is the possibility to define for any test function \I'(x,D) G L"^ the probability 
distribution p{E), with p{E)dE giving the probability to find the energy of the system in 
the interval [E,E + dE], such that the total probability integrated and/or summed over 
the whole energy range is unity. This is especially clear in the action-angle representation. 
The requirement is that, for every bounded function / and every test state \1/ G L^, the 
matrix element 



r+oa 

mXHsm = E / \Ml)\'f{ncr{I))dI 

= /(O) / \Ml)\'dI+ / \^n{I)\'f{ncr{I))dI. 



(43) 



is well defined (we have written the contribution of the isolated spectral point = as a 
distinct term). Now, '^n{I) are the Fourrier components of the test function \1/(J, </)). In 
original variables, their role is played by the integrals 

j ^{x,D)'^ENix,D)dxdD , (44) 

These integrals converge (though the normalization integrals for en{x, D) do not) and 
the weak singularity cx 1/ y/x — Xq does not hinder this convergence. 
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4.2 Including fermions 



Once the bosonic problem is resolved, it is not difficult to obtain the solution of the full 
problem ([6]). To be more precise, the Hamiltonian ([6]) as it is may appear problematic. For 
example, it does not look Hermitian — not invariant under the conjugation ijj ■ir^ ip , x ^ X 



Still, the spectrum of the Hamiltonian can be found by supersymmetry considerations. 
The states are classified by the value of the fermionic charge F, which can take values 
—1, 0, 1. The wave functions of the states in the sectors F = —1 and F = 1 involve the 
factor X and ipy correspondingly. The fermion part of the Hamiltonian does not act on 
such states and the solutions to the stationary Schrodinger equation in these sectors can 
be immediately written. 



^^=-\x,D;^,x) 
<f''=\x,D;^,x) 



(46) 



with written above in Eq.( |3711 . The eigenstates of the same energy in the sector F = 
can be obtained from the states (l46ll by the action of the supercharges Q,Q,T,T, which 
commute with the Hamiltonian. There are two such states. 



<ff=\x,D;ij,x) 



T-Q 



f + Q 



F=l 



^/2[N + V{x)] - iV\x)^IJX 



and 



F=0, 



x,D]'ip,x) 



Q + T 



D-^X-^ 
ox 



(47) 



^b{x,D). (48) 



The action of Q, T on and the action of T on \1>^=~^ give zero. 

And here we meet a certain difficulty. As far as the states (1471) are concerned, every- 
thing is fine, these states behave in the same way as involving a benign integrable 

are finite, and the prob- 



singularity oc 1/ ^Jx — Xq. The generalized Fourrier integrals 
ability for a test function \l/ from C2 to have an energy within a given interval {E, E + dE) 
is well defined. However, the functions fHS]) are more troublesome. They involve a singu- 
larity oc l/(a; — Xo)^^^, which is not integrable in the usual sense. 
We have two options now: 

1. To declare the state ( HHi) inadmissible by that reason. In that case, the spectrum 
would lose its supersymmetric form and would involve not quartets, but triplets of 
states of given energy, as illustrated in Fig|3l 



*Note, however, that it is invariant with respect to the involution 

^ X, X ^ -Va -0 ^ X, X ^ --0, 

supplemented by the usual complex conjugation of the bosonic variables. 



(45) 
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inadmissible 



inadmissible 



Q+T 



0>. T 



Q-T 




Q+T 



Q-T 



Q+T 



Q+T 



Q 



F = -l 



F=0 



F=l 



Figure 3: The triplet of "benign" states under the action of supercharges 



2. The second option is to include the states ( HHl) into some extended Hilbert space 
and to define the integrals (1^ with a kind of principal value prescription. Roughly 
speaking, we represent 



dx 



X 



3/2 



as 



- 2 / dx'^ix] 



d 1 



dx 



(49) 



(recall that the structure ~ l/x^^^ appeared in the first place when applying the 
differential operator to ^b)- We can integrate then by parts and write 



dx 



X 



with an integrable singularity. For mathematical precision, we can regularize 



1 



1 



1 



dy 



exp{-yVe^} 



do the integral in the R.H.S. of Eq. (H9|) at finite e and take the limit e ^ afterwards. 

The second option seems to us more reasonable. Once there is a consistent way to 
ascribe meaning to the states and the associated generalized Fourrier integrals (jUj) such 
that one of the important symmetries of the problem, the supersymmetry, is kept intact, 
it is preferrable to use it. In addition, the system of eigenfunctions involving only the 
states fH6l) and fH7|) . but not fHHj) does not seem to be complete!^ 



^We cannot formulate it as a rigourous mathematical statement, because we are dealing with contin- 
uum spectrum here, and the states (|46l47l48p do not represent a conventional basis in a conventional 
Hilbert space. The usual way to treat continuous spectrum systems is to put them in a finite box to 
make the motion finite and the spectrum discrete. The problem is that a nai've box regularization breaks 
supersymmetry. More study of this question is required. 
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As was mentioned before the Hamiltonian is not Hermitian. Still, its spectrum is 
real. That means that this Hamiltonian belongs to the class of so-called quasi-Hermitian 
or crypto- Hermitian Hamiltonians [12-14]. 



5 Models in the neighborhood. 

The Lagrangian ([2]) is the simplest nontrivial supersymmetric Lagrangian with higher 
derivatives. But there are many other such theories. In this section, we discuss two 
different natural modifications of (12]) 



5.1 Mixed theory. 

One obvious thing to do is add to ([2]) the standard kinetic term multiplied by some 
coefficient 7 and write 



L = J dOde 
The component expression for the Lagrangian is 



-{VX)^{VX) + -VXVX + V{X) 
2 dt 2 



(50) 



xb + DV'{x) + V"{x)xiIj + + 7 



x2 + D 



2 + ij^{i^x-i^x) 



The canonical Hamiltonian is convenient to express as 

h = h,-1f , 

where F is the operator of fermion charge and 



(51) 



(52) 



Ho = pP- DV'ix) -\{D^ + P^) + + 



(53) 



One can find out also the Nother supercharges Q, Q- Being expressed via canonical 
momenta, they are 



Q = ^[p + ^v'{x)]-(^x + li^){P-^D), 
Q = -x[p-^V'{x)] + (^i: + ^x){P + ^D) 



(54) 



Further, one can guess the existence of the following generalization for the second pair of 
the supercharges T, T, 



T 



^P[p-^V'{x)] + [x-li^) {P + iD) 



f = ^[p + rV\x)] + [i,-^-x]{P-iD). 



(55) 
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Introducing also the operators = xi^ and F_ = ifjx-, one can observe that the super- 
algebra of the set of the operators Hq, F, F_; Q,Q,T,T is closed. The nonvanishing 
(anti) commutators are 



[Q,Ho] = -|q, [Q,Ho] = |g, [T,Ho] = |t, [f,Ho] = -|f , 

[Q,F] = -Q, [Q,F] = Q, [T,F]=T, [f,F] = -f, 
[g,F_]=f, [Q,F+] = -T, [T,F^] = -Q, [f,F+] = Q, 
{Q,Q} = 2Ho-jF, {T,T} = 2Ho + jF, {Q,T} = 2jF+, {Q,T} = 2jF^. (56) 



One can make here a few remarks. 

• The operators F,F^,F- form the sl{2) subalgebra. One could have introduced the 
operators F± also in the case 7 = 0, but that was not necessary for closing the 
algebra. When 7 7^ 0, it is. Actually, Eq. (|56l) represents a well known simple 
superalgebra sZ(l,2) = osp{2,2) [15]. 

• The algebra fl56|) involves two conventional A/" = 1 subalgebras. They are realized by 
the subsets {Hq — (7/2)F; Q, Q) and {Hq + /2)F; T, T). Recall, however, that the 
operators Q, Q and T, T are not Hermitially conjugate to each other, which allows 
for the presence of the negative energies in the spectrum. 

• The algebra ( l56l) is a close relative of the unconventional so called weak supersym- 
metric algebra introduced in Ref. [16]. One can show that the latter is a semidirect 
sum of the algebra (!56|) with the Abelian 1-dimensional algebra {Y). 

• When 7-^0, the subalgebra of ( j56l) involving only the operators Hq = H, F; Q, T, Q, T 
coincides with the subalgebra of (|T5l) involving the same operators. The system with 

7 = involves an additional integral of motion A^, but, when 7 7^ 0, the Lagrangian 
is not invariant with respect to the transformations f|T6|) anymore and there seems 
to be no such integral. 

What is the dynamic of the mixed system ? Consider first the classical bosonic dy- 
namics. The Hamilton equations of motion are now 



The absence of the extra integral of motion A^ seems to make the system not integrable. 
That means that analytic solutions do not exist, but it is possible to study the solutions 

i°It is difficult to prove the absence of something. One can always say that the extra integral actually 
exists, but we simply have not found it. We performed, however, a numerical study which suggests that 
the system is not integrable. In particular, when 7 7^ 0, the parametric plot of the solution in the plane 
{x,P) does not represent a closed curve as it does for 7 = 0, but densely covers a certain region in the 
phase space. 



D 



P 
P 



X 



DV'\x) , 
V'{x) + 7D , 
P. 

p — 7P . 



(57) 
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numerically. Remarkably, it turns out that the trajectories are in this case in some sense 
more benign that for undeformed system. When 7 = 0, the function x{t) = XQcn[Qt,k] 
varies within a finite region, but the amplitude of the oscillations for D{t) grows linearly 
in time (see Eqs.fl22[ [25]) and Fig.l). For nonzero 7, it does not and the motion is finite. 
When 7 is small, the amplitude pulsates as is shown in FigJH The larger is 7, the less is the 
amplitude and the period of these pulsations. When 7 is large, the "carrying frequency" 
and amplitude fluctuate in an irregular way. 



D 




Figure 4: The function D(t) for a deformed system ( a; = 0, A = 1, 7 = .1). 

What can one say about the structure of the spectrum ? The finiteness of motion 
suggests that the spectrum might be discrete. Let us show that it is not discrete in the 
precise mathematical meaning of this word. A discrete spectrum involves only isolated 
eigenvalues [17]. We will prove, however, that the system under consideration involves an 
infinite number of eigenvalues in any finite energy interval. 

Let us consider first the sectors F = ±1 where the problem is equivalent to a purely 
bosonic problem with the Hamiltonian 

Hb = pP- DV'{x) - ^{D^ + . (58) 

Hb is unitary equivalent, up to a partial Fourier transform, to 

Kb = - {Pd, + V\x)dp) - 1{P^ - dl), (59) 

where dp = 

We can see that Kb has the same algebraic structure as the Fokker-Planck operator Ell 
but there is a big difference: the Fokker-Planck operator Hpp is 

Hl'p = yd^ + V\x)dy-^-{y'-dl) 

so Hpp is not Hermitean but the Hermitean part of Hpp is negative (7 > 0) or positive 
(7 < 0) definite. Here Kp is Hermitean and is unbounded both from above and from 

^^We thank B. HelfFer for this remark. 
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below, and we have lost the "hypoelliptic character" of the Fokker-Planck operator (see 
[18] for more details). 

Nevertheless, for harmonic potentials V{x) = we can compute the spectrum 

of Hb explicitely. To diagonalize the (classical) Hamiltonian 0581) . we introduce first 
canonically conjugated holomorphic variables 

p — iu'^x ^ p + iuj'^x . P — iD _ P + iD 
a = a = ^ ; b = ^— , b = ^— , (60) 

such that the Hamiltonian acquires the form 

Hb = uj{ba* + ab*) - 

Then we rotate 



Ci \ _ / COS0 sm0 \ / a 
C2 J I — sin cos (j) ) \ b 



(61) 



with tan(20) =2uj/'-j = t. The Hamiltonian is reduced to 

Hb = uJiCicl - UJ2C2C2 , (62) 



where 



The spectrum 



tOi,2 = l[VlT^Tl 



ujjk = (^uji + J - (^^2 + ^ + const, j, k eN, (63) 

of the quantum counterpart of this Hamiltonian is pure point ( a constant shift in ( l63l) 
is inserted for generity. It takes into account the ordering ambiguity). Moreover, this 
spectrum is dense if uji/uj2 is not rational, which is true for generic r. 

It is worth noting that the system (158!) with quadratic V{x) is equivalent to Pais- 
Uhlenbeck oscillator — a higher-derivative dynamic system described by the Lagrangian [3] 

L = I {if - {ul + + ululq) (64) 

The Hamiltonian can be found using Ostrogradsky formalism [19]. If quantizing it and 
evaluating the spectrum, we obtain exactly the same result as in Eq. (1631) [5,20]. 

Remark: For more general potential V , it is not immediately clear that Hb is es- 
sentially self- adjoint. This can be proved if V"{x) is uniformly bound on M. But we do 
not know how to obtain rigourous results about the spectrum of Hb for non quadratic 
potentials. If the coupling 7 is imaginary, we are in the Fokker-Planck situation and, 
if some confining conditions on the potential V are satisfied, it is known then that the 
spectrum of Kb is pure point with bounded from below eigenvalues. Moreover, in this 
case the resolvent of Kb is compact [18]. 
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For real positive 7 and more general potentials V a semi-heuristic reasoning revealing 
the nature of the spectrum can be presented. Consider the quantity Z[f] = TT{f{HB)}, 
where f{u) is any positive definite function dying at m = ±00 fast enough. For example, 
one can take f{u) = exp{— m^/o"^}. In semiclassical approximation, we can evaluate it as 

Z[f]^l'^^^mt(.,D,p,P)], (65) 

where H'^ is the Weyl symbol of the quantum Hamiltonian [in our case, it is given directly 
by Eq. fl58l) ]. The corrections to this formula [their existence can be understood by noting 
that the Weyl symbol of /(-ffs) does not coincide with /(if^^)] can also be evaluated [21]. 
When the function f{u) is smooth enough (for f{u) = exp{— M^/a^} the condition is 
(7^1), the corrections are small. Doing in Eq. flB^ the integral over dDdP, we obtain 

with 

giu) = / f{w)dw . 



When X and/or p and hence u are large, the integrand in fl66l) is a constant, and the 
integral diverges. On the other hand, assuming that the spectrum is pure point, one may 
write 

Z[f] = (67) 

n 

The infinite value of this sum for any function / including the functions that die at infinity 
very fast means the presence of an infinite number of states in a finite energy range. That 
means that the spectrum should have accumulation points so that the eigenvalues of the 
Hamiltonian are not well separated from each other. It is conceivable that in our case each 
point of the spectrum is an accumulation point and the spectrum represents a countable 
subset of M that is dense everywhere. This is what happens for simpler models with the 
Hamiltonian like 

^ = ^ + ra!- V + D^) (68) 
27 27 2^ ^ ^ ^ 

(the semiclassical value of Z[f] for the Hamiltonian H is the same as for Hb). 

Another possibility is that the spectrum of Hb is truly continuous with not normaliz- 
able wave functions. Based on the mentioned above fact that the classical motion of our 
system is finite, we find this option less probable. But only a future study will allow one 
to obtain a definite answer to this question. 

Let us briefiy discuss the dynamics of the full supersymmetric system. As was ex- 
plained above, it involves two pairs of complex supercharges. However, the supercharges 
T, T do not commute with the Hamiltonian H = H_ = Hq — {'j/2)F, but only with the 
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operator = Ho + {•y /2)F. On the other hand, Q and Q do not commute with H+. The 
fact that the maximal conventional supersymmetry subalgebra involves in this case only 
one pair of complex supercharges implies that the spectrum of H consists of degenerate 
doublets rather than quartets (as was the case for 7 = 0). The same concerns H^. 

The doublet structure of the spectrum is a feature which distinguishes the system 
under consideration from the system considered in [16]. The algebra of the latter was 
similar to fl56l) . but involved an extra bosonic charge Y. That allowed for the existence 
of an operator that commutes with all supercharges. The spectrum of this operator (it 
is natural to call it Hamiltonian) beyond the ground state and the first excited state is 
4-fold degenerate. 



5.2 More derivatives. 

As a final example, consider a somewhat more complicated action 



S 



dtdOde 



-iVX) iVX) + ViX) 



The corresponding component Lagrangian is 

The bosonic equations of motion are 

x^^^+DV"{x) : 
D-V'{x) 

In the simplest quadratic case 

ViX) = ^ 



+ ^^piJ + DV'{x) + V"{x)iIji{j . 








(69) 



(70) 



(71) 



(72) 



the equations (ITTll are linear and can readily be solved. Their characteristic eigenvalues 
are 



Ai,: 



X 



3,4,5,6 




(73) 



We see that, besides oscillating solutions, there are also solutions with exponentially 
growing amplitude. Strictly speaking, the Hamiltonian is still Hermitian in a certain 
sense due to the fact that there is no collapse: it takes an infinite time to reach infinity 
and a unitary evolution operator can be defined at all times. However, Hermiticity is lost 
as soon as one switches on interactions. We solved numerically the equations of motion 
(I7T!) for V{x) oc ±x^ and V{x) oc ±x^ and found out that the solutions collapse reaching 
a singularity at finite time. 



-'^^The commutation relation [T,H] = '^T guarantees that, if ^' is the eigenstate of H, is also 

an eigenstate, but with a different eigenvalue. Thus, only the supercharges Q, Q are effective as far as 
degeneracy is concerned. 

^^modulo a possible remedy in the spirit of [13], see the discussion below. 
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6 Discussion 



Probably, the main lesson to be learned from the analysis of different higher-derivative 
quantum mechanical models in this paper is that the ghosts (negative energy states and 
the Hamiltonians without bottom) do not always lead to violation of unitarity, but one 
should worry about it only in the case when the collapsing classical trajectories exist. The 
analysis performed in Refs. [1,6] displays that sometimes even in this case the quantum 
problem is (or can be, if defining the Hilbert space with a care) well defined, but for the 
model ([2]) where there is no collapse, quantum evolution is unitary in spite of the absence 
of the ground state. 

In addition to this, we found a bunch of rather unusual phenomena. 

1. The classical trajectories of the system ([2]) do not collapse, but exhibit oscillatory 
behavior with linearly rising amplitude. For the modified model ( l50l) . the amplitude 
does not grow and the trajectories are finite. 

2. The model ([2]) is exactly soluble due to the presence of an extra integral of motion. 
For the quartic potential, the solutions of the classical and quantum problems are 
expressed analytically via elliptic functions. 

3. Besides Nother supercharges, the systems ([2]) and flSU]) involve an extra pair of 
supercharges. Nother supercharges, the additional supercharges, the Hamiltonian 
and certain extra operators form a modified supersymmetry algebra. For different 
systems, these modified superalgebras are also different. 

4. The system ([2]) has continuous spectrum. The spectrum of the model (|5Ui) is prob- 
ably not continuous, but involves only normalizable discrete spectrum states, with 
a countable set of eigenvalues densely covering M. This conjecture needs to be 
confirmed. 

5. The spectrum of the Hamiltonian (jH]) was found to be real even though this Hamil- 
tonian is naively not Hermitian in the fermion sector. That means that it belongs 
to the class of so called quasi-Hermitian, alias crypto-Hermitian Hamiltonians hav- 
ing attracted recently a considerable interest [12-14]. A salient feature of such 
systems is the possibility to define a modified norm in Hilbert space such that the 
Hamiltonian is manifestly Hermitian with respect to this norm [27]. It would be 
interesting to find out whether such modified norm exists also in our case. 

The central question posed in Refs. [1,2], whether benign higher- dimensional higher- 
derivative supersymmetric field theories exist, is still left unresolved. The best currently 
known candidate for this role is superconformal at the classical level renormalizable gauge 

^^They are unusual for conventional systems with positive definite kinetic term, but maybe not so 
unusual for the system involving ghosts. Unfortunately, the latter were never seriously studied before. 

^^In supersymmetric context, crypto-Hermitian systems were first discussed in [22]. It was found 
recently that the Hamiltonians describing so called nonanticommutative supersymmetric theories [23] 
are crypto-Hermitian [24,25]. Non-Hermitian supersymmetric a models, having certain kinship to our 
models, were considered recently [26]. 
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theory in six dimensions constructed in Ref. [28]. Unfortunately, it is not free of difficulties 
and it is not clear at the moment whether they can be resolved or not. The simplest such 
model constructed in [28] is probably not viable because of chiral anomaly. However, the 
anomaly is cancelled in a theory involving besides gauge supermultiplet also a matter 
hypermultiplet [29]. On the other hand, the latter scale- invariant theory is not fully 
conformal even at the classical level and has an infinite number of propagating fields [30] . 

In addition, in contrast to the models discussed in the present paper, the models [28,30] 
do involve collapsing classical trajectories reaching infinity in finite time. This is due to 
the presence of the cubic term ~ D"^ in the Lagrangian. (Fields D are the highest 
components of the vector M = 1 QD supermultiplet of canonical dimension 2. They are 
auxiliary for the standard quadratic in derivatives theory, but become dynamical when 
extra derivatives are added.) The toy models analyzed above provide no care for such 
a theory. The latter still might exist in some sense. Indeed, if the theory with complex 
cubic potential V{x) oc ix^ involves a benign crypto- Hermitian Hamiltonian [13], the 
theory with real cubic potential oc may also acquire sense if considering it on the 
complex x-plane 

It is not inconceivable that the higher derivative theory of Ref. [28] can also be treated 
along these lines. On the other hand, if we are interested not just in quantum mechanical 
systems, but in field theory, we would like to have not only Hermitian Hamiltonian, but 
also unitary S'-matrix. A discouraging news [25] is that, for the field theories involving 
the cubic term in the potential, the unitarity of S'-matrix is impossible to preserve, if 
choosing the asymptotic states in a conventional way. It is not clear, however, whether 
we have to require the existence of conventional asymptotic states and the conventional 
unitary S - matrix for the fundamental theory in the higher-dimensional bulk. Maybe, it 
would suffice to have a unitary finite time evolution operator with a certain (complicated) 
choice for the Hilbert space metric ? 

Further studies of all these questions are necessary. 

We are indebted to E. Ivanov, B. Helffer, V. Kac, M. Kroyter, N. Nekrasov, and S. 
Theisen for illuminating discussions and correspondence. 
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